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Foreward 


This  report  is  Chapter  XI  of  the  twelve  in  a  forthcoming  research 
monograph  on  the  mathematical  theory  of  laminar  combustion-  Chapter  I-IV 
originally  appeared  as  Technical  Reports  Nos.  77,  80,  82  &  85;  these  were 
later  extensively  revised  and  then  issued  as  Technical  Summary  Reports 
No's  1803,  1818,  1819  8t  1888  of  the  Mathematics  Research  Center,  University 
of  Wisconsin-Madison.  References  to  I-IV  mean  the  MRC  reports. 


Contents 

1.  Scope 

2.  Slowly  Varying  Perturbations  of  the  Plane  Wave 

3.  Buoyancy  and  Curvature 

k.  Stability  of  Near-Eauidiffusional  Plane  Flames 

5.  Cellular  Flames 

6.  Hydrodynamic  Effects 
7-  Curved  Cellular  Flames 

References 


page 

1 

3 

10 

lit 

20 

26 

30 

37 


Figures  1-9 


The  findings  of  this  report  are  not  to  he  construed 
as  an  official  Department  of  the  Army  position  unles 
so  designated  by  other  authorized  documents. 


Chapter  XI 
Stability 


1»  Scope 

Flame  stability  is  a  subject  of  enormous  breadth,  which  could  easily 
fill  a  monograph  by  itself  albeit  a  less  analytical  one  than  ours;  both 
the  mechanisms  and  manifestations  of  instability  assume  many  forms.  It 
has  been  studied  at  least  since  the  observations  of  Higgins  (1777)  on 
singing  flames,  a  phenomenon  involving  interaction  of  the  kind  described 
in  Chapter  V  between  an  acoustic  field  and  an  oscillating  flame.  Such 
interactions  were  laboratory  curiosities  for  many  years,  but  lately  they 
have  assumed  technological  importance  in  the  development  of  rocket  motors 
and  large  furnaces.  This  type  of  instability  is  distinguished  by  being  well 
understood  (Chu,  ) ,  at  least  qualitatively.  Other  types  of  instability 
are  in  general  poorly  understood  and,  in  some  cases,  even  in  doubt,  in  the 
sense  that  instability  is  suspected  of  playing  a  role  but  no  certain 
evidence  is  yet  available.  The  following  examples,  though  by  no  means  con¬ 
stituting  a  complete  list,  convey  some  idea  of  the  richness  of  instability 
phenomena  in  combustion. 

The  stability  of  burner  flames  depends  on  an  appropriate  interaction 
between  flame  and  surroundings,  in  particular  a  flux  of  heat  from  the  flame 
to  the  burner  rim.  The  role  of  this  flux  in  anchoring  the  flame  and  prevent¬ 
ing  blow  off  has  already  been  mentioned  in  Chapter  IX. 

Propagation  limits  are  important  questions  for  premixed  flames,  steady 
sustained  combustion  being  possible  only  for  a  certain  range  of  the  fuel  to 
oxidant  ratio.  If  a  mixture  is  too  rich  or  too  lean  it  will  not  burn,  and 
instability  may  well  play  a  role  in  the  matter.  A  possible  related 
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phenomenon  is  quenching  by  heat  loss  or  non-uniformities  (Chapter  X),  but 
whether  this  is  a  question  of  existence  of  stability  (or  some  combination  of 
the  two),  is  not  yet  known. 

Flow  is  turbulent  in  the  right  circumstances  and  indeed  turbulent  com¬ 
bustion  is  of  paramount  importance  in  many  technological  applications.  How 
burning  influences  the  transition  from  laminar  to  turbulent  flow  is  not  well 
understood  at  present.  Instability  may  also  play  a  role  in  another  transi¬ 
tion,  that  from  deflagration  to  detonation  (DDT).  A  flame  propagating  through 
a  mixture  will,  under  appropriate  stimulus,  accelerate  and  change  from  being 
an  isobaric  deflagration  wave  into  a  detonation  wave,  with  shock-like 
structure. 

Instability  is  invariably  invoked  to  eliminate  branches  of  multiple 
responses  for  being  physically  unrealizable.  Several  examples  are  provided 
by  the  burning  of  a  condensate  (Chapter  IV) ,  and  others  by  fuel-  and  mono¬ 
propellant-drop  combustion  (Chapters  VI  and  VII). 

Finally,  a  striking  manifestation  of  instability  is  provided  by  cellular 
flames,  which  are  discussed  in  detail  by  Markstein  (1964).  Under  certain  con¬ 
ditions  a  nominally  plane  wave  displays  cells,  i.e.  an  unsteady,  quasi-periodic, 
two-dimensional,  transverse  pattern  (Fig.  1).  On  average  the  cells  have  a 
characteristic  dimension,  but  large  cells  grow  and  divide  while  small  cells 
shrink  and  vanish.  The  corresponding  phenomenon  in  highly  curved  burner 
flames  leads  to  beautiful  polyhedral  tips  (Fig.  2),  ->hich  sometimes  are  quite 
Steady  and  sometimes  rotate  but  never  exhibit  the  unsteadiness  of  plane 
cellular  flames.  Cells  have  also  been  observed  on  a  spherical  flame  spread¬ 
ing  out  from  an  ignition  source  (Istratov  &  Librovich  1969). 


All  these  examples  apply  to  premixed  flames  and  some  to  diffusion  flames, 
which  introduce  their  own  peculiar  questions  associated  with  the  separation 
of  the  reactants.  However,  we  shall  be  concerned  only  with  the  stability 
of  the  freely  propagating  premixed  flame,  and  cellular  flames  will  play  a 
central  role  in  the  development.  Such  a  narrow  focus  is  dictated  by  the 
state  of  the  art,  analytically  speaking:  certainly  it  is  the  only  area  that 
has  been  explored  by  activation-energy  asymptotics  to  any  siginificant  extent, 
a  big  attraction  being  the  simple  closed  form  of  the  unperturbed  steady  state. 
Even  then  formidable  difficulties  remain;  that  they  have  ben  overcome  is  due 
to  the  tenacity  with  which  Sivashinsky,  the  principal  architect  of  the  theory' 
we  shall  describe,  has  pursued  the  problem. 

In  our  opinion,  the  material  in  the  present  chapter,  answering  as  it 
does  questions  that  have  long  frustrated  the  combustion  scientist,  is  one 
of  the  triumphs  of  activation-energy  asymptotics. 

2.  Slowly  Varying  Perturbations  of  the  Plane  Wave 

The  motion  of  a  slowly  varying  flame  is  governed  by  the  basic  equation 
(VIII. 32),  where  v^  is  to  be  determined  from  a  coupled  hydrodynamic  analysis. 
There  the  flame  appears  as  a  discontinuity  across  which  the  jump  conditions 
(VIII. 1,2, 3)  hold  and  on  either  side  of  which  Euler  »s  equations  (VIII. 34)  are 
valid. 

If  the  gas  is  quiescent  for  ahead  of  the  flame  there  is  a  solution 

-1  “  0  »  Mn  =  1  »  K  =  0  »  (D 

namely  the  steady  plane  flame.  Small  perturbations  of  this  stationary  state, 
to  be  denoted  by  primes,  then  satisfy 
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where 


M!  =  p/v'-F')  , 

n  1  nl  t 


T*Pl/b 


and  we  have  used  V  =.  1/p^.  Without  loss  of  accuracy,  the  £,  r),  5-frame 
may  be  fixed  at  an  arbitrary  point  of  the  undisturbed  plane  wave.  The 
velocity  perturbation  must  be  determined  by  the  hydrodynamic  analysis 
mentioned  above,  a  step  to  which  we  shall  return  shortly. 

Only  for  one-dimensional  disturbance  of  the  flame  surface  can  the 
velocity  perturbations  be  neglected,  as  will  become  clear  in  the  sequel. 
Nevertheless,  doing  so  in  general  leads  to  a  correct  description  of  one  of 
the  disturbance  modes  as  k+0  or  |g  [•><».  (The  latter  corresponds  to 
small  heat  release,  when  the  unit  of  temperature  becomes  small  and  hence  T* 
large,  or  to  a  Lewis  number  close  to  one.)  Sivashinsky  (1977a)  then  speaks 
of  the  dif fusional- thermal  effect,  since  these  are  the  processes  left.  With 


v'  =  0  and  hence  M'  =  -p,F'  , 
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equation  (2)  admi^ts  spatially  periodic  solutions 


F'  =  PjA  exp  (ax  +  ikn) 


provided 


(k^-p^a2-28p1a)A 


jjCt  =  -6±  /e2  +  k2  , 


0  i.e. 


(6) 


m 
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where  no  loss  of  generality  results  from  considering  disturbances  Independent 
of  ^ . 

For  one-dimensional  disturbances  k  =  0  the  two  roots  are 

=  0  ,  -20  ,  (7) 

and  the  flame  is  unstable  if  and  only  if  0  is  negative,  i.e.  the  Lewis 
number  is  greater  than  1.  This  result  was  uncovered  in  our  earlier  discussion 
of  plane  flames  (Sec. HI. 5).  On  the  other  hand,  for  k  ^  0  there  is  a 
positive  root  whatever  the  value  of  £  so  that  instability  always  occurs. 

The  conclusion  .is  found  to  be  correct  when  velocity  perturbations  are  taken 
into  account,  although  the  smaller  (in  absolute  value)  of  the  roots  (6),  on 

which  it  is  based,  is  spurious,  i.e.  not  associated  with  any  of  the  actual 

^fu 

modes.  On  the  other  hand.^larger  root  does  correspond  to  an  actual  mode  in 
the  limit  |  0  |  ~>  °°  or  k  -*•  0,  though  only  to  leading  order.  Thus  there  is 
a  diffusional-thermal  mode  but  it  is  not  necessarily  unstable.  We  turn  now 
to  the  hydrodynamic  problem. 

The  jump  conditions  (VIII. 1,2, 3)  show  that  the  state  behind  the  plane 
flame  is 


p2  =  pl“Yl  ’  u2  =  Y1  »  v2  =  0  <8> 

while  the  perturbations  on  the  two  sides  satisfy 

p^  =  p^-2(o-l)(u’-F;)  ,  =  au’-(o-l)F’  ,  v^  =  v^F’  . 

(9) 

Here  u  and  v  are  the  velocity  components  in  the  and  p-  direction 


and 
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o  =  12/Tj^  =  1  +  Y1/T1 


(10) 


When  account  is  taken  of  the  convective  accelerations  due  to  the  velocities 
T^(=  l/p^)  and  T2  in  front  of  and  behind  the  discontinuity,  Euler's 
equations  (VIII. 34)  are  seen  to  have  the  solutions 


[p'.u'.v']  =  [-  (p  j^a+|  k|),|k|  ,ik]P  ^exp(ctt+|  k  |  £+ikq)  for  £  <  0  , 

[p'.u'.v']  =  [p2a-|k|  ,  |  k  |  ,-ik]P2exp(aT-|k[£+ikri) 

+[0,o|k|  ,-ip^asgnk]Sexp(ax-p2a£+ikri)  for  £  >  0 


The  P-terms  form  potential  fields  for  the  velocities  while  the  S-terms,  which 
are  only  solenoidal,  take  account  of  (convected)  vorticity  behind  the  flame. 
These  solutions  are  appropriate  for  a  flame  surface  (5),  and  substituting 
into  the  jump  conditions  (9)  yields 


o[p^a+(2a-l)  |  k|  ]P^  +  (p^a-crj  kj  -  2a(a-l)p^aA  =  0  , 

a|k|P1-  ]kj P2  -  cr| k]  S  -  (a-Dp^A  =  0  , 

I  kj  Px  +  |k|P2  +  p1aS-  (a-l)IkjA  =  0  . 


If  the  flame  displacement  A  were  specified,  the  equations  would  determine 
the  hydrodynamic  fields  on  the  two  sides  of  the  surface.  In  fact  A  is 
related  to  for  slowly  varying  flames,  as  the  basic  equation  (2)  shows. 

Landau  (1944)  and,  independently,  Darrieus  (1946)  supposed,  in  the 
absence  of  a  theory,  that  the  flame  speed  is  unaffected  by  the  perturbations; 

'tfu 

which,  as  we  shall  see,  is  true  in^theory  of  slowly  varying  flames  for 
k  0  or  |$|  -*■  ®.  The  response  of  the  flame  is  thereby  eliminated  and  it 
is  appropriate  to  speak  of  the  hydrodynamic  effect.  Setting  the  flame-speed 


y 


perturbation  -  F^  equal  to  zero  gives 

|k|P1-p1aA  =  0  ,  (13) 

and  then  the  homogeneous  system  for  P^,  P2,  S  and  A  has  a  non-trivial 
solution  if  and  only  if 

2  2  0 

PjCt  =  a  |  k  |  or  (o+l)p^a  +2o  |  k|  p  ^a-a  (a-1)  k“  =  0 

i.e.  pa  =  a(-l ±  /1+a-l /c) , k  / (0+ 1)  .  (14) 

The  first  root  gives  P  =  A  =  0  so  that  the  stability  of  the  flame  is  not 
involved,  although  disturbances  behind  the  flame  are  amplified.  However,  A 
is  non-zero  for^ the  second  root,  which  is  alvavs  positive  in  view  of  c  >  1, 
and  we  conclude  that  the  flame  is  unstable.  Note  that  the  result  is  quite 
independent  of  the  theory  of  slowly  varying  flames  (even  though  £,  n,  C 
were  used):  the  hydrodynamic  effect  is  always  destabilizing  though,  of  course, 
it  may  only  be  present  when  the  variations  are  slow. 

We  now  come  to  the  stability  of  the  plane  wave  to  slowly  varying  dis¬ 
turbances,  including  both  the  dif fusional-thermal  and  hydrodynamic  effects. 

The  basic  equation  (2)  provides 

(k2  +  Pla|ki  +2$|k|)P1+ (k2-p^a2- 2Bp1a)A  -  0  (15) 

as  fourth  equation  of  the  homogeneous  system,  in  place  of  (13),  so  that  the 
condition  for  a  non-trivial  solution  becomes 


0 


(16) 
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where 


a  =  o+l  ,  a. 
o  1 


a3 


(o+l) ( ( k|+2B)  ,  a2  =  o|k|[-(3o-l)|k|+4B] 

-ok2[ (3o-l) |kj+2(a-l)B]  , 


* 

J 


(17) 


if  p^a  =  a { k |  is  again  discarded. 

From  their  behavior  as  k  •+  0  or  |  B  |  -*■  »  it  is  clear  that  one  of  the 
three  modes  may  be  called  dif f usional-thermal  and  the  other  two  hydrodynamic. 
In  either  case  the  cubic  becomes 


(p^a+26)  [  (o+l)p^a2  +  2ojk|  p^a  -  o(o-.l)k2]  =  0  (18) 

so  that  the  larger  of  the  roots  (6)  and  the  two  relevant  roots  (14)  are 
obtained.  Correspondingly  the  basic  equation  (15)  reduces  to 

(p2a2 +  2Bpja)A  =  0  .  (19) 

That  the  condition  (19)  agrees  with  the  dif f usional-thermal  condition  (6) 
only  in  the  limit  k  -*■  0  or  [  g  |  -*•  00  itself  shows  that  the  smaller  root 
(6)  is  spurious.  It  predicts  p^a  =  0(k  /2B)  for  k  small  or  |8j  large, 
whereas  the  actual  remaining  modes  have  p^a  =  0(k). 

The  cubic  shows  that  steady  corrugations  (a  =  0)  with  wavenumber 

|k|  =  -2(c-l)B/(3o-l)  (20) 

are  possible  when  S  is  negative  (£  >  1) .  The  result  was  first  obtained 


by  Sivashinsky  (1973)  but  Ls  of  historical  interest  only  since  for  that  value 
of  k  (and  indeed  all  values)  we  shall  see  that  there  is  at  least  one  root 
of  the  cubic  with  positive  real  part,  implying  instability. 


V 
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A  result  of  the  type  (16),  for  a  more  complicated  model,  was  first 
obtained  by  Eckhaus  (1961)  from  a  thin-flamc  analysis  which  did  not  use 
formal  activation-energy  asymptotics.  Instead  he  adopted  a  formula  of  Mallard 
and  LeChatelier  (see  Emmons,  1958*  which  leads  to  the  perturbations  of  flame 
temperature  and  speed  being  proportional,  a  relation  concistent  with  the 
asymptotic  result  (VIII. 30)  after  linearization.  The  first  complete  deriva¬ 
tion  was  given,  somewhat  heuristically ,  by  Sivashinsky  (1977b)  and,  rigorously, 
by  Buckmaster  (1977) 

Necessary  conditions  for  every  root  of  the  cubic  (16)  to  have  a  negative 
real  part  are  a^,  a^  >  0.  Stability  would  therefore  require 

-M  <  26  <  -(^)|ki 

which,  for  o  >  1  as  here,  is  a  contradiction.  We  are  forced  to  conclude 
that  plane  flames  are  unstable. 

The  instability  occurs  whatever  the  wavenumber  k,  but  it  would  be  a 
mistake  to  infer  that  all  steady,  slowly  varying  flames  are  unstable.  Such 
an  extrapolation  implies  taking  the  short-wave  limit  k  -*■  °°  (so  as  to  screen 
out  boundedness  and  curvature)  and  hence  violating  the  assumption  of  slow 
variations  on  which  the  conclusion  of  instability  for  unbounded  plane  flames 
was  based.  In  any  event  the  consequences  of  linear  instability  are  not 
necessarily  catastrophic;  the  result  may  be  a  nonlinear,  unsteady  structure 
superimposed  on  the  steady  state  without  obliterating  its  essential  character. 

Even  the  instability  of  unbounded  plane  flames  to  long-wave  disturbances 
cannot  be  considered  definitive,  since  there  may  be  stabilizing  influences 
that  have  been  ignored  so  far.  Two  such  influences  are  gravity  and  curvature. 


which  we  now  examine. 


The  influence  that  gravity  has  on  flames  by  virtue  of  the  induced 
buoyancy  forces  is  clear  from  a  commonplace  observation:  an  inverted  candle 
does  not  behave  in  a  desirable  fashion.  A  less  obvious  manifestation  is  the 
dependence  of  the  vertical  propagation  limit  (i.e.  the  dilution  beyond 
which  the  mixture  is  not  inflammable)  on  whether  the  propagation  is  upwards 
or  downwards.  Instability  should  play  a  role  in  such  limits  since  for  up¬ 
ward  propagation  the  light  burnt  gas  underlies  the  heavy  unburnt  gas,  so  that 
Taylor-type  instability  can  be  anticipated.  (The  characterization  is  not 
entirely  accurate  since  the  flame  is  not  a  material  surface.) 

A  flame  of  characteristic  dimension  1cm.  and  speed  45  cm. /sec.  has 
approximately  equal  gravity  and  dynamic  heads.  These  quantities  are  typical, 
and  faster  or  smaller  flames  are  correspondingly  less  influenced  by  gravity. 
Moreover,  since  flame  thickness  is  typically  much  less  than  a  millimeter  (Sec. 
II. 4),  gravity  can  only  be  important  in  the  hydrodynamic  field  of  a  nominally 
plane  flame.  Support  for  that  conclusion  comes  from  the  derivation  of  the 
basic  equation  (VIII. 32)  of  slowly  varying  flames,  which  does  not  use  momentum 

balance  (the  only  governing  equation  influenced  by  gravity,  see  Sec.  1.2). 

3  5 

The  unit  of  acceleration  is  cpM  /ocA,  which  is  typically  2  *  10  cm. /sec 

-3 

The  dimensionless  gravity  force  is  then  5  *  10  and  we  may  reasonably  con¬ 
sider  it  0(1/9).  The  momentum  equation  (VIII. 34)  now  gains  a  term  -p(^,0) 
on  its  right-hand  side,  where  g  is  positive  for  downward  propagation,  so 
that  terms 


-PjgS  and  -p2g£ 


(21) 


are  added  to  the  pressures  in  front  of  and  behind  the  flame,  respectively. 
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The  perturbations  and  p^  in  equation  (9)  are  now  the  sums  of  terms 

coming  from  equation  (11),  as  before,  and  terms 

-pjgF’  =  -p^g  A  exp (aT+ikn)  ,  -p,,g  F*  =  -p^  ^  exp(ax+ikr\)  , 

so  that  the  first  of  conditions  (12)  is  replaced  by 

a[p1a+(2o-l)  |  k|  ]P^+ (p^-o  |k.|  )P2  -  (o-l)  (2ap1a-p^g)A  =  0  .  (22) 

Clearly  there  is  no  limiting  effect  on  the  dif fusional-thermal  mode  but 
there  is  on  the  hydrodynamic  modes,  the  roots  (14)  becoming 

P-,0  =  o  [-1  ±  /l+(o-l/o)  (1-p^g/a  |  k  1 )  1 1  k |  /  (o+l)  .  (23) 

Since  g  is  negative  for  upward  propagation,  we  immediately  see  that  gravity 
is  destabilizing  when  the  hot  gas  underlies  the  cold.  When  the  reverse  is 
true,  disturbances  for  which 

M  <  p^g/a  (24) 

are  stable.  For  a  downwardly  propagating  flame,  modes  of  sufficiently  long 
wavelength  are  therefore  stabilized. 

That  these  conclusions  do  indeed  hold  is  seen  from  the  cubic  (16)  with 

terms 


(o-l)p^g| k|  (p jCt  +  |k[  +  26) 

added  to  account  for  gravity.  In  either  of  the  limits  k  -*■  0  or  8  -*■  00 


the  cubic  becomes 
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(p^ct+26)  [  (o+l)p^a2  +  2a|k|p1a  -  a(a-l)k2  +  (c-l)p2gjk!  ]  =  0  (25) 


so  that  the  diC fusional-thermal  mode  is  unaffected  but  the  modifications 
(23)  of  the  hydrodynamic  modes  are  obtained.  For  the  latter,  the  stabilizing 
effect  is  most  clearly  seen  for  long  waves,  when 


P 


o 

a-1 


|k|  +0(k3/2) 


(26) 


We  conclude  that  for  £>  i  gravity  cannot  completely  stabilize  even  long 
waves,  since  the  relevant  root  (7)  is  a  creature  of  purely  one-dimensional 
disturbances  and  these  are  unaffected  by  the  hydrodynamics;  but  for  £  <  1 
it  can. 

As  k  -*■  <=°  the  cubic  reduces  to 

(Pla+|k|)[(o+l)p2a2-a(3o-l)k2]  =  0  ,  (27) 

which  is  independent  of  g:  short  waves  are  unaffected  by  gravity,  as 
expected.  Note  how  instability  persists  in  the  limit:  there  is  always  a 
positive  root,  though  the  other  two  are  negative. 

We  now  examine  the  effect  of  curvature,  but  only  on  the  limit  of  the 
dif fusional-thermal  mode.  A  full  treatment  is  quite  complicated,  much  more 
so  than  for  a  plane  flame,  and  in  any  event  gravity  stabilizes  the  other  modes 
(at  least  for  long  wavelengths).  Attention  will  be  confined  to  a  plane  source 
or  sink  flow 

u  =  c/r  ,  v  =  0  (28) 

where  9r  is  radial  distance  and  u,  v  are  polar  components.  The  undisturbed 
flame  is  then  circular,  with  curvature 
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k  *  -1/p1c  (29) 

determined  by  the  requirement  that  the  mass  flux  is  1  at  it. 

Such  a  flame  corresponds  to  the  solution 

vpl  =  0  ,  Mn  =  1  ,  V  =  0  (30) 

of  the  basic  equation  (VIII. 32),  and  perturbations  of  this  steady  state 
satisfy 


ft’  -  V,-v|,  -  kV'  +—  M*  =  0  . 

*1 


n  -J_  -pi  p,  n 


(31) 


Here  v'^  is  due  to  the  displacement  of  the  flame  in  the  source  or  sink 
flow,  disturbances  of  which  have  been  ignored.  It  is  sufficient  to  consider 
plane  disturbances,  for  which  the  flame  may  be  written 


-  *  P-Jcl  +  F'  (<)),t) 


(32) 


in  polar  coordinates  r,  4> .  Then 


3  2 


-  -  *f;  ,  m;  -  TtFVojc+p^;)  ,  ,  (33> 


according  as  c  is  positive  or  negative,  and  the  behavior  of  modes 


F'  =  PjA  exp(aT+in<fi) 


(34) 


is  determined  by 


p^a2  +  2(8+l/p1c)p1a  +  (2B-n2/p1c) /p^c 


0 


(35) 
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as 


/  "  ■>  ~ 

=  -(6+1/Pj^c)  i  *  (ii+l/p^c)*' -  (2B-n_/p1c)/o1c  .  (36) 

Clearly  the.  roots  are  always  real. 

To  discuss  the  result  we  introduce  the  wavenumber 

k  =  u/pjjcl  (37) 

and  note  that  the  formula  (6)  is  recovered  in  the  limit  ]cj  ®  with  k 
fixed.  Since  the  smaller  root  was  then  found  to  be  spurious,  we  shall 
discard  the  corresponding  root  here.  The  larger  root  is  negative  if  and  only 
if  6  >  -1/p^c,  so  that  curvature  is  a  stabilizing  influence  for  the  source 
flow  but  destabilizing  for  the  sink  flow.  The  conclusion  must  however  be 
considered  tentative  until  incorporation  of  the  hydrodynamics  shews  that  we 
are  in  fact  dealing  with  the  limit  6  “  or  k  -*■  0. 

The  requirement  that  c  be  positive  corresponds  to  a  flame  that  is 
concave  towards  the  fresh  mixture,  as  for  a  closed  burner  flame.  Such 
flames  were  found  for  £  <  1  (Sec.  IX. 3)  and  our  theory  predicts  stability 
of  the  dif fusional-thermal  mode.  On  the  other  hand,  the  open  flames  found 
for  £  >  1  correspond  to  c  <  0  and  instability  is  inferred.  Extrapolations 
of  this  kind  are,  however,  little  more  than  guesses. 

4.  Stability  of  Near-Equidif fusional  Plane  Flames. 

The  conclusion  that  modes  of  short  wavelength  will  always  be  unstable 
according  to  the  theory  of  slowly  varying  flames,  while  not  a  valid  one, 
nevertheless  stands  in  sharp  contrast  to  experimental  results.  We  therefore 


turn  to  near-equidLf fusional  flames  for  clarification,  and  here  compromise 
is  necessary.  On  the  scale  of  the  flame  thickness  slowly  varying  flames  are 
almost  plane  and  almost  steady;  these  characteristics  allow  the  fluid 
mechanics,  in  particular  variable  density,  to  be  incorporated  without  difficulty. 
By  contrast,  the  continuity  and  momentum  equations  do  not  simplify  when  the 
flame  is  near-equidif fusional,  so  that  perturbations  of  even  the  plane  flame 
lead  to  non-constant  coefficients.  For  that  reason  most  discussion  of  near- 
equidif  fusional  stability  has  relied  on  the  constant-density  approximation, 
as  embodied  in  equations  (VIII. 62,  63),  and  with  p  single  exception  of  Sec.  6 
that  is  the  framework  adopted  here. 

For  flames  in  a  quiescent  gas,  the  governing  equations  reduce  to 

3(T,h)/3t  =  V2 (T ,h+XT)  ;  (38) 

to  these  are  added  the  jump  conditions  (VIII. 64).  The  steady  plane  wave  cor¬ 
responds  to  the  solution 


T  =  Ti  +  Yje”  ,  h  *  -VJ^ne11  for  n  <  0  , 

T  =  ,  h  »  0  for  n  >  0 

where 


(39) 


n  =>  x+ 1  , 


(40) 


the  flame  sheet  being  located  at  x  =  -t  [cf.  equations  (IX.29,30)].  Its 
stability  to  small  perturbations  is  our  first  concern. 

The  class  of  disturbances  is  limited:  0(1)  perturbations  of  T  are 
admitted  provided  they  are  balanced  by  0(1)  perturbations  of  Y  which 
keep  the  sum  of  order  9  \  It  follows  that  a  finding  of  stability  (in 
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contrast  to  instability)  is  not  definitive.  Nevertheless  the  conclusions 
are  very  convincing  and  it  seems  unlikely  that  an  analysis  dealing  with  a 
larger  class  of  disturbances  would  alter  the  picture  substantially. 

Such  perturbations  satisfy 

ad’.h^/at  =  v2(t'  ,h'+XT')  (41) 

on  both  sides  of  the  disturbed  flame  sheet 

n  =  F’(y,t)  .  (42) 

The  appropriate  solution  is 

T'  =  0  for  n  >  0  (43) 


behind  the  flame  sheet.  Continuity  of  T  and  h,  plus  the  jump  conditions 
on  their  derivatives,  yield 


T'  =  -YjF'  ,  [Vi’  ]  =  -XYjF'  , 


3T'/3n  =  -Yj^F’ +Y1h'/2T2  ,  [3h'/8n]  =  XY^’  /2T*  -  2XYXF' 


for  n  =  0 


where  the  right-hand  sides  of  the  last  two  equations  are  to  be  evaluated  for 
n  =  0+,  i.e.  on  the  hot  side  of  the  flame  sheet. 

Variations  in  the  z-direction  may  be  ignored  without  loss  in  generality 
of  the  conclusions.  The  solution  corresponding  to 

F’  -  A  exp  (at  +  iky)  (45) 


is 


(T * , h * )  =  (B,C+X (i<2-k2)Bn/ (1-2k^) )exp(at+iky-h<^n)  for  n  <  0  , 

(46) 

h’  =  D  exp  (at+Lky-h^n)  for  n  >  0  , 

where 

2  =  yfli  ^ L+4ct+4k2 ]  ,  -tt  <  arg(l+4a+4k“)  <  -  .  (47) 

To  be  admissible  the  solution  must  vanish  far  ahead  of  the  flame  sheet  and 
be  bounded  far  behind,  i,e.  Re  must  be  positive  and  Re  <2  non-positive.. 
For  Re  a  >  0  (i.e.  the  unstable  modes  of  interest  and  their  neutral  limits) 

these  requirements  are  met  automatically,  so  we  need  not  concern  ourselves 
with  them  further.  Equations  (44)  ncH' yield  a  homogeneous  system  for  4,  B, 

C,  D  which  has  a  non-trivial  solution  only  if 

(l-K1)(l-2ic1)2  +  X[(l-<1)2-k2]  =  0,  X  =  YjX/212  ,  (48)  „ 

a  result  first  obtained  by  Sivashinsky  (1977a). 

This  condition  determines  all  those  modes,  i.e.  values  of  a,  correspond¬ 
ing  to  any  given  X  and  k.  For  Re  a  >  0  we  can  be  sure  the  mode  is 
admissible,  the  boundaries  of  regions  of  instability  Re  a  >  0  being 
curves  Re  a  =  0.  There  are  just  three  such  curves,  namely 

the  X-axis:  k  *  0  ,  (49) 

the  parabola  P:  4k2  =  -(X+l)  ,  (50) 

the  branch  B  of  the  curve  (l+12k2)X2  -  4(l+8k2)X  - 
-  8<l+8k2)2  =  0  for  which  X  >  16/3  . 


(51) 


For  the  first  two  Im  a  =  0  also. 

The  regions  of  instability  can  now  be  determined  from  explicit  results 
for  small  k,  namely 

ry  »  /•  rt  n 

ct  =  - (X+l)k  +  X  (X-3)k  +0(k  )  ,  [X  -4X-8+x/x(X-8)]/32  +  0(k  ) 


where  both  square  roots  are  implied  for  0  <  X  <  8.  Clearly  there  is  a  mode 
with  Re  a  >  0  for 


X  <  -1  or  >  2(l+/3) 


respectively,  but  not  otherwise,  so  that  the  stability  regions  are  as  marked 
in  Fig.  3.  Note  how  misleading  results  would  be  obtained  by  considering  only 
one- dimensional  disturbances:  for  k  =  0  the  first  of  the  values  (52)  vanishes 
and  no  instability  is  predicted  for  X  <  -1.  The  resulting  stability  condition 
X  <  2(l+/3)  can  be  considered  a  refinement  of  jf  <  1  for  slowly  varying 
disturbances  [Sec.  2,  cf.  equation  (7)]. 

For  each  value  of  k  there  is  a  band  of  Lewis  numbers,  always  including 
-1  <  X  <  16/3  and  hence  =  1,  for  which  the  flame  is  stable.  The  conclusion 
is  in  accord  with  that  for  slowly  varying  disturbances,  where  instability  is 
found  for  all  X  /  1  (Sec.  2)  and,  moreover,  provides  a  possible  explanation 
of  why  stable  flames  are  observed,  a  possibility  that  makes  near-equidif fusional 
flames  of  great  practical  importance.  To  be  sure,  in  the  limit  8  «  the 

stability  band  in  £  is  vanishingly  small,  but  in  practice  9  is  only 
moderately  large.  Since  Lewis  numbers  are  invariably  close  to  one,  it  is  con¬ 
ceivable  that  all  flames  adequately  described  by  activation-energy  asymptotics 
may  be  considered  near  equi-dif fusional. 


y 
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On  the  boundary  B  the  neutral  stability  modes  are  oscillatory,  i.e. 

In  a  ^  0.  The  boundary  is  relevant  to  mixtures  whose  deficient  component 
is  heavy  and  therefore  has  relatively  small  diffusivity  (e.g.  for  lean 
hydrocarbon  or  rich  hydrogen  flames  in  air),  but  it  has  received  no  attention 
to  date. 

On  the  parabolic  boundary  and  on  the  X-axis  the  neutral  modes  are  non- 
oscillatory,  i.e.  Im  a  =  0.  (In  fact,  all  modes  inside  the  parabola  are 
non-oscillatory  too.)  For  any  X  <  -1  and  wavenumber  in  the  range 


0  <  k  <  /-(X+D/2 


there  is  a  mode  that  will  grow.  Moreover,  since  RecC  vanishes  at  the  ends 
of  the  range  there  must  be  an  intermediate  value,  namely 


.4+151-91  +  (4-6X)/l-3X.,l/2 
i  1  108X  J 


for  which  the  growth  is  a  maximum.  Although  the  instability  does  not  favor 
one  wavelength  to  the  exclusion  of  all  others,  an  arbitrary  disturbance  can 
be  expe.cted  to  grow  in  such  a  way  that  the  length  1/k^  plays  a  significant 
role.  As  we  noted  in  Sec.  1  it  is  characteristic  of  cellular  flames,  which 
are  usually  observed  when  the  deficient-component  is  light  (X  negative) ,  for 
many  modes  to  contribute  to  the  structure,  but  nevertheless  an  average  dimension 
can  be  assigned  to  the  cells. 

That  the  stationary  mode  corresponding  to  the  parabola  has  relevance  to 
cellular  flames  is  supported  by  the  behavior  of  the  flame-temperature  perturba¬ 


tion 


20  - 


2T^(1-k^)A  exp  iky  =  T^(A+4k^  -  1)A  exp  i(ky-hi)(k  =  /-(X+l)/2) 


corresponding  to  the  displacement 

F*  =  A  exp  iky  ,  (57) 

with  which  it  is  exactly  out  of  phase.  Viewed  from  the  burnt  mixture,  the 
crests  are  colder  than  the  troughs.  We  would  therefore  expect  the  crests  to 
be  the  least  luminous  part,  which  is  in  fact  a  well-known  characteristic 
of  cellular  flames. 


5.  Cellular  Flames 

The  linear  analysis  of  the  previous  section  predicts  that  infinitesimal 
disturbances  will  grow  without  bound  when  X  is  less  than  -1,  but  there  is 
str^ong  evidence  that  the  growth  actually  leads  to  cellular  flames.  In  that 
case,  nonlinearities  must  limit  the  growth;  the  nature  of  that  interaction 
will  now  be  investigated.  Attention  will  be  focussed  on  the  neighborhood 
of  X  =  -1  by  setting 


-X  I  e 
oT 


2T*/Y! 


where  e  >  0  is  vanishingly  small.  Only  the  top  sign  on  z  is  of  interest 

here,  corresponding  to  the  linearly  unstable  waves  in  that  neighborhood 

(i.e.  the  long  ones),  but  the  bottom  sign  will  be  used  later.  Such  a  focus 

can  be  treated  by  a  perturbation  analysis  in  e.  The  stability  boundary  (50) 

1/2 

show-  that  these  waves  have  a  wave-number  k  =  0(e  )  at  most,  while  the 

2 

result  (52)  shows  that  they  grow  no  faster  than  a  =  0(e  ).  Accordingly  the 


slow  variables 
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(59) 


are  appropriate  for  their  description. 

For  fixed  c,  an  infinitesimal  disturbance  will  grow  according  to  the 
linear  analysis  until  its  size  makes  some  nonlinear  term  comparable  to  the 
linear  terms.  Its  amplitude  is  then  0(c)  so  that,  following  Sjrtishinsky 
(1977b),  we  take  the  flame  sheet  as 


x  =  -t  +  eF(n,r) 


(60) 


Nonlinearity  has  certainly  intruded  by  then,  as  is  seen  from  the  resulting 
flame  speed 


W  =  1-c3(Ft+  F3/2)+0(e6) 


(61) 


but  it  could  conceivably  intrude  elsewhere  for  amplitude  o(e);  the  consistency 
of  the  analysis  stemming  from  the  assumption  (601  shows  that  it  does  not. 

If  now  the  coordinate 

1/7  7 

n  =  x+t~eF(e  y,c  t)  (62) 

is  introduced  in  place  of  x,  so  that  n  =  0  is  the  flame  sheet,  then  derivatives 
in  equations  (38)  become 


_3_  _  _3_  _3_  3/2  J_  1/2  _3_  _3_ 

3x  3n  ’  3y  6  n  3n  E  3n  *  3t 


while  the  normal  derivative  in  the  jump  conditions  is 


(l-c3F  )-3-+c2  £ - 


u+!  e3f2)  -  c2f  -L 

1  +  2  e  V  3n  E  Fn  3n 
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Perturbation  expansions  are  now  introduced  for  T,  h  and  F,  for 
example 


T 


2  3 

T  +  cT  +  e  T.  +  e  T,  +  . .  , 
o  1  2  3 


(63) 


where  the  coefficients  are  functions  of  n,  n,  t.  A  sequence  of  problems 
then  results,  whose  solutions  are  subject  to  the  requirements  that 

=  •••  =0  for  n  >  0,  that  conditions  as  n  ■+  -<*>  are  undisturbed 
and  that  exponential  growth  as  n  -+  +°°  is  disallowed  (though  algebraic  growth 
cannot  be  prevented).  At  the  fourth  problem,  for  and  h^,  the  nonlinear 

equation  (75)  for  Fq  is  obtained  as  a  solvability  condition. 

The  first  problem  is 


[TJ 


i<V-V2V3"2  ■  0 

[h  ]  =  0  ,  X  [3T  /3n] 

o  o  o 


for  all 


=  [ 3 h  /3n] 
o 


-X  Y  exp (h  Hi.) 
o  i  o  * 


where 


L=- 


2  2 
3  /3n 


3/3n 


(65) 


As  expected,  its  solution  is  the  steady  plane  wave 


T 

o 


+v 


C X  Y,ne 
°  1 

hQ  =  J  for  n  £  0  .  (66) 

l  0 


no  confusion  will  result  from  the  double  use  of 


T^,  since  we  shall  show 


immediately  that  the  expansion  coefficient  is  identically  zero.  In  fact  the 
second  problem  is  easily  seen  to  be 
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JCCTX)  =  0  , 

/(hL)  -  X^T^D 

"2  =  ft*' 

for  n  i  0  ,  \ 

V<(67) 

[Txl  =  [h1l 

=  0  ,  [3T  /3n] 

■  -"lA0  • 

[31^/30]  =  -(hj^  i  Y^)  at 

when  the  Independence 

of  T  and  h 

o  0 

from  n,  t 

is  used;  so  that 

1  (o  , 

r«i-* 

ni‘(o 

for  n  i  0 

(68) 

is  obtained  without  using  the  jump  conditio  on  3h^/3n,  which  is  then 
automatically  satisfied.  The  information  obtained  so  far  enables  the  third 
problem  to  be  written 


r<v  ■  ;v  ■  l<h  2>-vV3"2  • 


[T9]  =  [h7]  =  o  ,  [3T9/3n]  =  -h2/Xo  ,  [3h?/3n]  =  -h2 


Again,  the  solution 


fY,¥  ne 

1  onn 


T2  = 


f\  Y.F  (l+n2)en 
1  o  1  onn 


h2  =< 


(70) 


X  Y.F 
o  1  onn 


is  determined  without  using  the  jump  condition  on  Dh^/Sn,  which  is  then 
satisfied  automatically. 

Finally  we  reach  the  fourth  problem,  which  can  now  be  written 


-  (Y^a.n+a^e  r  2 

jC(T,)  =  )  .  -C(h3)  -  Xq32T3/3 


n2  =fY l<a3n 

KVl 


2  .  .  n 

+a^n+a^)e 


for  n 


[T  1  -  [h3l  =  0  ,  [DT3/3n]  =  Y1Fon/2-h3/XQ  ,  Oh^n]  =  -h3*Vonn+X 


Note  that  must  be  allowed  to  increase  algebraically.  The  jump  condition 

on  Sh^/Sn  is  not  satisfied  automatically  but  requires 

F  +4F  i  F  /X  +F2  /2  =  0  .  (75) 

or  onnnn  onn  o  on 

This  result  is  due  to  Sivashinsky  (1977b) ,  who  generalized  it  to 

F  +4V4F  ±  V2F/>.  +  (VF  )2/2  =  0  (76) 

OT  O  O  '  O 

foundations  in  both  the  y-  and  z-  direccions. 

The  equation  has  an  interesting  implication  for  the  flame  speed.  While 
the  expression  (61)  is  purely  kinematical,  its  approximation 
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W 


1+  e3(4F  ±  F  /X  ) 

onnnn  onn  o 


(77) 


is  a  profound  consequence  of  the  dif fusional- thermal  processes.  Markstein 

(1964),  in  an  attempt  to  modify  Landau's  conclusion  that  all  flames  are  hydro- 

dynamically  unstable,  assumed  that  the  flame-speed  perturbation  is  proportional 

to  curvature  and  so  incorporated  the  term  in  F  but  not  the  other. 

1  onn 

The  linearized  form  of  equation  (75)  has  solutions  proportional  to 
exp(ar/e^+ikn/c1^)  if 

a  =  ek^/X  -  4k^  ,  (78) 

o 


which  coincides  with  the  first  of  approximations  (52)  in  the  neighborhood  of 

X  *  -X  ,  k  =  0  •  This  result  is  graphed  in  Fig.  4.  The  maximum  growth  rate 
o 

occurs  for 


k.  =  /T78X-  ,  (79) 

i  o 

which  is  an  approximation  to  the  general  formula  (55).  The  connection  of 
this  maximum  with  the  characteristic  cell  size  has  already  been  noted. 

Existing  discussion  of  the  nonlinear  equation  (75)  is  limited  to  numerical 
computations  (ifeivashinsky//&(>Michelson),  1977) .  Integration  of  the  initial- 
value  problem  with  periodic  boundary  conditions  (using  a  period  large  compared 
to  k^)  leads  to  results  like  those  in  Fig.  5  when  the  plus  sign  is  taken 
in  the  equation.  (For  the  minus  sign  there  is  decay  to  zero,  corresponding 
to  linear  stability.)  The  nonlinearity  prevents  unlimited  growth  of  the 
disturbance  and  the  resulting  structure,  being  quasiperiodic  with  characteristic 
dimension  k^,  highly  suggestive  of  cellular  flames. 


6.  Hydrodvnamic  Effects 

The  analysis  of  Sec.  5,  leading  as  it  does  to  a  balance  of  small  effects, 
provides  a  framework  within  which  several  generalizations  can  be  discussed  so 
long  as  the  additional  effects  are  appropriately  small.  Here  we  shall  con¬ 
sider  weak  interaction  of  the  flame  with  the  hydrodynamic  field.  A  rigorous 
derivation  of  the  governing  equation  replacing  (75)  requires  analysis  similar 
to  that  in  Sec.  5.  A  more  transparent  derivation  can  be  made  by  plausible 
arguments . 

The  constant-density  model  is  studied  because  it  is  simple  and  yet 

retains  some  of  the  key  physical  ingredients  that  govern  flame  behavior. 

Nevertheless  it  is  a  rational  asymptotic  limit  (a  -*•  1)  ,  albeit  one  of  limited 

practical  interest  since  density  changes  across  flames  are  rarely  small.  When 

o  is  close  to  1,  the  hydrodynamic  effects  can  be  incorporated  in  a  rational 

fashion  as  perturbations  just  strong  enough  to  influence  the  flame  behavior 

_1/? 

described  in  Sec.  5.  Since  we  are  only  concerned  with  wavelengths  0(t  “)  , 

the  field  outside  the  flame  sheet  has  a  dual  structure  similar  to  that  of 

slowly  varying  flames  (Sec.  VIII. 3):  hydrodynamic  disturbances  described  on 

-1/2 

the  scale  n  =  0(e  )  must  be  matched  with  those  in  the  diffusion  zone, 

where  n  =  0(1) . 

The  flame  speed  (61)  is  now  augmented  by 

(u;  -  e3/2F  v;)[l  +  0(t3)l 
i  n  i 

where  uj,  are  the  components  of  the  velocity  distui banco  immediately 

ahead  of  the  flame  sheet.  If  this  Is  to  be  aporturbat 1 on  just  comparable  to 

3 

the  existing  one,  then  u^  must  be  0(c  );  in  which  case  v^  (being  of  the 

3 

same  order)  makes  a  smaller  contribution  and  V,T  is,  to  order  e  ,  simply 


augmented  by  u^. 
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The  problem  now  is  Co  express  u'  in  terms  of  F.  which  involves 

0  ® 

treating  the  flame  sheet  in  the  same  way  as  the  hydrodynamic  discontinuity 
was  in  Sec.  2.  Indeed,  the  weakness  of  the  velocity  perturbation  ensures 
that  Euler’s  equations  are  satisfied  even  in  (he  diffusion  zone,  so  that  we 
may  read  off  the  result 


PL  =  (o-l)A/2  (80) 

of  solving  the  jump  conditions  (12)  when  k  and  a  are  appropriately  small 

and  0  is  close  to  1.  Equation  (80)  relates  u^  and  F0  when  they  are 

proportional  to  exp(at+iky)  and  hence  the  Fourier  transforms  of  general  u| 

and  F  (since  a  is  not  involved).  It  follows  that 
o 


dk 


1 


,  n  3/2  F(n,T) 
(g-De  /  on _  j- 

2»Pl  -L  ^  d  ’ 


(81) 


where 

a  =  l  +  0(e3/2)  (82) 

if  is  to  be  0(e3). 

So  far  we  have  been  concerned  with  modifying  the  kinematical  result  (61) 
for  the  flame  speed;  now  we  come  to  the  result  (77)  of  the  dif fusional-thermal 
processes,  which  is  found  to  be  unchanged.  Certainly  it  is  plausible  that 
such  weak  density  changes  do  not  affect  these  processes.  That  being  the  case, 
the  generalization 
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and  Fig.  6  shows  graphs  of  these  curves.  As  expected  from  the  last  section, 
for  y  =  0  the  lower  sign  gives  stable  modes  only,  but  the  upper  sign  involves 
unstable  inodes.  For  y  >  0,  i.e.  when  the  hydrodynamic  field  is  disturbed, 
both  signs  give  unstable  modes.  The  possibility  of  stabilizing  the  inodes 
again  by  gravity  can  also  be  seen.  Gravity  changes  the  velocity  by  an  amount 
proportional  to  the  displacement  of  the  flame  sheet  (cf.  Sec.  3),  which  cor¬ 
responds  to  adding  a  constant  to  the  righ-hand  side  of  equation  (85).  If 
the  constant  is  negative  (i.e.  gravity  points  from  burnt  to  unburnt  gas)  and 
sufficiently  large,  then  the  curves  in  Fig.  6  are  translated  downwards  into 
the  quadrant  a  <  0. 

Sivashinsky  &  Michelson  (1977)  have  treated  the  full  equation  (83)  also, 
in  particular  for  the  pjifynj  sign  when  the  hydrodynamics  is  an  extra  destabiliz¬ 
ing  effect.  Results  similar  to  Fig.  5  are  obtained,  except  that  now  the 
large  wrinkles  themselves  contain  several  small  wrinkle-.,  '.,'e  shall  consider 
just  two  limits  of  the  equation,  both  corresponding  to  significantly  larger 
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V 


(89) 


30  - 


The  equation  expresses  constancy  of  the  flame  speed,  as  Landau  assumed  in 
his  stability  analysis. 

A  continuous  periodic  solution  for  F  ,  obtained  numerically,  is  show 

L4tu.  V/y  '.1%^ 

in  Fig.  7/T  (A  solution  for  any  other  period  can  be  obtained  from  it  by 
j’U,-’.  jgcufc  £ /uijw  Vi' 

scaling  nj[ )  The  <i  is  continuities  in  Fq^  at  ±tt  are  smoothed  out  by  the 
neglected  *-orir)-term>  whose  inclusion  perturbs  the  flame  speed.  We  are 
therefore  led  to  the  suggestion  that,  under  more  general  circumstances,  the 
destabilizing  effect  of  hydrodynamics  first  discovered  by  Landau  does  not 
result  in  chaotic  structure  but  just  a  wrinkling  of  the  nominal  flame  con¬ 
figuration.  Such  fine  wrinkles  are  sometimes  observed  on  actual  flames 
(Markstein,  1964). 

7.  Curved  Cellular  Flames. 

The  cellular  pattern  in  Fig.  1  is  formed  by  a  flame  located  in  a  large 
(10  cm.  diameter)  vertical  tube  and,  consistent  with  the  results  of  Sec.  5, 
is  highly  unsteady.  By  contrast,  cellular  flames  stabilized  on  burners  are 
usually  steady,  which  suggests  a  classical  bifurcation  phenomenon  similar  to 
Taylor  cells  in  cylindrical  Couette  flow.  At  a  critical  value  of  some  para¬ 
meter  disturbances  of  a  single  finite  wavenumber  become  unstable,  and  beyond 
that  value  there  is  an  alternative  (stable)  steady  state.  Here  we  shall 
establish  curvature  as  such  a  parameter,  leadingjto  a  steady  cellular  pattern. 
As  in  Sec.  3,  attention  will  be  confined  to  the  plane  flow 

u  =  c/r  ,  v  =  0  (90) 

with  c  now  positive  only  (source  flow);  r  is  measured  on  the  scale  of 
the  flame  thickness.  If  disturbance  of  the  hydrodynamic  field  is  again 
neglected  then  equations  (38)  govern  provided  the  replacement 
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3  +  J_+ c  _3_ 

3t  3t  r  8r 

'is  made.  These  possess  a  steady  solution  satisfying  the  jump  conditions 
(VIII. 64)  in  which  the  flame  sheet  is  located  at 

r  =  c  (91) 


and 


T  =  T^  +  Y^(r/c)C  ,  h  =  XY^(rC/cC  ^)in(c/r)  for  r  <  c  ,  ^ 

\  (92) 

T  =  ,  h  =  0  for  r  >  c  .  J 

Such  is  the  undisturbed  flame. 

The  result  (91)  implies  that  the  dimensionless  burning  rate  is  1,  i.e. 

the  dimensional  burning  rate  is  the  same  as  that  of  an  adiabatic  plane  f.lame. 

There  is  an  apparent  contradiction  with  previous  sections,  where  we  saw  a 

change  in  burning  rate  with  curvature,  which  is  resolved  on  recognizing 

change  in  curvature  as  the  determining  factor.  The  essential  process  is 

transverse  diffusion,  which  is  absent  from  the  cylindrical  flame  because  of 

symmetry,  and  its  presence  is  detected  by  changes  in  curvature.  The  simplified 

treatment  to  be  presented  next  will  ensure  that  only  disturbances  in 

curvature  influence  the  burning  velocity. 

The  objective  is  to  admit  curvature  disturbances  just  strong  enough  to 

2 

modify  equation  (75).  Since  the  F^-term  corresponds  to  curvature  0(e  ), 
we  are  therefore  led  to  introduce 


2 

e  c 


(93) 


w 


(99) 


x 
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where  v  is  here  the  velocity  in  the  y-direction.  However,  the  radial 

flow  is  locally  just  a  weak  simple  strain  superimposed  on  a  uniform  flow, 

so  that  the  result  (99)  is  recovered  by  letting  the  strain  rate  become  small 

in  Sec.  X.3.  It  is  no  accident  then  that  X  =  -X  plays  a  role  both  in 

o 

the  effect  of  straining  flows  on  a  flame  and  in  linear  stability  analysis. 

A  modification  of  the  kinematical  result  (61)  arises  because  the  flame 

sheet  moves  in  a  non-uniform  velocity  field  (without  disturbing  the  field). 

3  3/2 

The  accompanying  changes  (-e  F^/R,  e  n/R)  in  fluid  velocity  at  the  sheet 
lead  to  the  new  approximation 


W  =  1  -  £3 [F  +F2  / 2  +  R~ 1  ( F  -  nF  )] 

or  on  o  on 

=  l-e3(f  +f2/2R2  +  f/R)  . 

T  0) 


(100) 


Comparing  the  results  (98)  and  (100)  establishes 


f  +4R~4f 

X  OJOIUOI 


±  R  2f  /X  +R~1f  +  R-2f2/2 

OMj)  O  U) 


0 


(101) 


as  the  fundamental  equation,  governing  the  radial  perturbation  f  in  the 
position  of  the  flame  sheet.  The  term  R  ^f  arises  from  the  perturbation 
in  radial  velocity  at  the  flame,  so  that  for  nominally  spherical  flames 
(where  the  velocity  is  proportional  to  r  )  the  term  is  replaced  by  2R  f 
(Sivashinsky ,  1978).  We  shall  pursue  the  stability  of  the  cylindrical 
flame;  the  spherical  flame  requires  only  minor  changes,  including  generaliza¬ 
tion  of  the  m-derivatives  (which  is  also  needed  here  for  variations  in  the 
z-direction) . 

For  linear  stability  we  set  f  proportional  to  exp(ax+ivuj)  and  neglect 
the  quadratic  term  to  obtain 
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a  =  ±R  2\)2/X  -4R-4y4-R1 
o 


(102) 


where,  although  it  is  strictly  a  discrete  parameter  (being  e 


3/2 


times  an 


integer),  v  effectively  takes  all  real  values.  (If  variations  in  the 

2 

z-direction  are  admitted,  with  wave  number  k,  then  v  is  replaced  by 

2  22  * y/ 

v  +R  k  and  the  conclusions  are  modified  in  an  obvious  way.)  As  R  -*■  °=  ‘Urt^k 


R_1V  =  k 


(103) 


held  fixed  the  result  (78)  is  recovered;  and  finite  R  clearly  has  a 
stabilizing  effect.  For  X  >  -X  curvature  only  argments  the  existing 
stability,  but  for  X  <  -X  it  can  overcome  the  instability  at  sufficiently 
small  wavenumbers  k  (see  Fig.  8) .  For  all  R  the  maximum  growth  rate  is 
attained  for  the  value 


1/8X 


(104) 


cf .  equation  (79)  ,  and  it  is  zero  when  R  has  the  value 


R  =  16X 

c  o 


(105) 


for  R  <  R  there  is  stability,  but  for  R  >  R  there  is  a  band  of  unstable 
—  c  c 

waves . 

2 

The  stability  regions  in  the  plane  of  /\Jk|  and  R/Xq  are  shown  in 
Fig.  9.  As  R  increases  through  Rc  a  small  band  of  unstable  wavenumbers 
appears,  suggesting  a  classical  bifurcation.  The  corresponding  Landau 
equation  (cf.  Matkowsky,  1970)  is  obtained  in  the  usual  way  be  setting 


e2  =  (R-R  )/R  ,  f  =  e(f  +  ef  +  e“f„ -i- . . . )  ,  t  =  x*/c2 

c  c  o  1  l 


(106) 
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with 


f  »  e  A  cosv  in  for  t*  =  0  (107) 

c 

To  leading  order  we  find  that 

f  =  A  (t*)cos(v  oi)  ,  A  (0)  =  A  (108) 

o  o  c  o 

and  then  that  secular  terms  appear  in  f2  unless 

A'  =  A  (9~A2A2)/144X2  .  (109) 

o  o  o  o  o 

It  follows  that  the  amplitude  asymptotes  to  the  constant  value 

A  =  ±3/A  according  as  A  $  0  (110) 

o  o 

We  conclude  that  the  cylindrical  flame  sheet  is  stable  for  R  <  but 

unstable  for  R  >  R  ;  and  that  for  small  positive  values  of  R  -  R  there 
c  c 

1/2 

is  a  stable  cellular  configuration,  with  amplitide  proportional  to  (R-Rc)  , 
to  which  the  originally  cylindrical  flame  sheet  tends  after  disturbance. 

A  characteristic  feature  of  stationary  cellular  flames  (including 
polyhedral  flames)  is  that  the  crests,  i.e.  the  portions  of  the  flame  that 
are  convex  towards  the  burnt  gas,  are  sharper  than  the  troughs.  Here  the 
variations  in  curvature  are  given  by 

o  /  0  *x 

f  *  e{-32A  [A  cos(v  u)+eA1cos(v  u)+eB1sin(v  w)]-(64X  /9)eA  cos(2v  u>)}+0(e  ) 
(iXl)  o  o  c  l  c  1  c  o  o  c 

(111) 

where  Aq  has  its  ultimate  value  (110)  as  do  the  integration  functions 

A, (t*)  and  B. (t*)  in  f,.  For  A„  >  0  the  crests  lie  at  os  -  2mr/v 
1  1  X  o  c 
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and  the  troughs  at  w  =  (2n+l)  it /\j ,  where  n  is  an  integer.  Clearly  the 

2 

0(c  )  harmonic  increases  the  magnitude  of  the  curvature  at  the  crests  and 
decreases  it  at  the  troughs,  consistent  with  observations. 
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